Transitions in quantum networks 
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We consider transitions in quantum networks analogous to those in the two-dimensional Ising 
model. We show that for a network of active components the transition is between the quantum and 
the classical behaviour of the network, and the critical amplification coincides with the fundamental 
quantum cloning limit. 
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Precise control over single quantum systems is essential 
in testing and harnessing of quantum mechanics. This 
has become possible with the advances in laser cooling 
and trapping techniques and manipulation of optical ele- 
ments in the one-photon level. The availability of single 
quantum systems has fed the interest in quantum net- 
works: A quantum computer Q is a network of individ- 
ual quantum systems, where any two of the nodes can 
interact with each other. Most quantum private commu- 
nication schemes Q are networks of two or three nodes. 
In addition to these information processing and commu- 
nication related applications, networks of optical com- 
ponents H and avoided crossings in multilevel systems 
|^,^) have been considered in order to study higher di- 
mensional quantum interference effects. 

According to statistical physics, a set of probabilisti- 
cally behaving individual systems can exhibit critical be- 
haviour when connected. In this paper we consider the 
question whether transition phenomena exist in networks 
of systems which behave probabilistically not because of 
finite temperature but due to their quantum nature; tran- 
sitions are known to exist in Ising quantum chain models. 
We define a model of a quantum network which carries 
in its structure a formal analogy to the two-dimensional 
Ising model. Such networks can be experimentally re- 
alized by various active (energy-consuming) or passive 
(energy-preserving) components. It is found that transi- 
tions do take place and we are able to give them a clear 
physical interpretation. For active systems the transition 
is between quantum and classical, for passive systems be- 
tween diabatic and adiabatic behaviour of the network. 
The transition phenomenon is clearly reflected in observ- 
able quantities. 

The Ising model describes a set of two-state systems 
which interact with their nearest neighbours; a quantum 
analogy of such a setup can be experimentally realized 
in various ways, as will be explained below. Fig.l shows 
schematically a 2-D quantum network with nearest neigh- 
bour interactions. To define the building blocks of this 
network we now take a closer look at the Ising model. 

The two-state systems in the Ising model, let us say 
spins, are on a 2-D lattice of the size N x M. Since 
only nearest neighbour interactions are taken into ac- 
count, the total energy of the system can be expressed 



using the energy E(fi) of one column (with the spin 
configuration /x) and the energy E(/i, //) between two 
columns. Let Sk = ±1 denote the values of individual 
spins and e be the absolute value of the energy of a spin- 
spin interaction. Then the energies can be written as 
E{n) = -eX)fc=i SfcSfc+i and E(/j,,fj,') = -eJ2k=i s ks'k- 




FIG. 1. A quantum network where the nodes are connected 
to their nearest neighbours. The boxes A and B denote the 



transformations performed at the nodes; when they are de- 
fined by Eq.(l) the network has a relation to the Ising model. 

The partition function Q(T) can be expressed in a sim- 
ple form by defining a 2 N x 2 N matrix V whose matrix 
elements are the thermal weight factors corresponding 
to a particular spin configuration of two neighbouring 
columns (n\V\n') = e -P\ E (w')+ E M\ . With this nota- 
tion Q{T) = TrV M = Eali( A ") M - The eigenvalues A Q 
thus determine the thermodynamics of the system. As 
was shown by Onsager and Kaufmann jfj), the 2 N x 2 N 
matrix V is a spinor representation of a set of plane ro- 
tations in 2JV x 27V dimensional space. The eigenvalues 
of V are uniquely determined by the eigenvalues of the 
corresponding plane rotation matrix P, which is 



P 



A 













A 



A 



B22 





B12 







B21 



B 



B 



where 



1 



A 



cosh 9 i sinh 6 
-i sinh. 6 cosh 9 

cothfl 



B 



cosh</> 
-i sinh( 



i sinh ( 
coshrf 



cosh( 



(1) 



and 4> = 2e/kT. 

The form of the matrix P suggests immediately a 
quantum-network analogue. The matrices A and B can 
be interpreted to describe unitary evolution of a two-state 
or two-mode system. The matrix P is then the evolution 
operator over a period in the network of Fig. 1 . The in- 
puts of the network are mixed pairwise according to the 
transformation A, and then the pairs are let to interact 
with the neighbouring ones by applying the shifted set of 
operations B. By repeating this M times, a 2N x 2M 
dimensional network can be constructed jjj . Since P con- 
tains all the physical information of the Ising model, e.g. 
the phase transitions, we may expect analogous phenom- 
ena in the quantum network described by P. Note that 
our aim is not to consider a set of quantized spin-systems 
at a finite temperature, like in the context of NMR quan- 
tum computing S. Instead, we borrow from the Ising 
model the abstract structure describing classical statisti- 
cal behaviour, and ask what kind of quantum behaviour 
it could describe. 

The physical realizations of the quantum network P 
can be divided into two groups. When the angle 4> 
is real, A and B are SU(l,l)-type matrices describ- 
ing energy-consuming (active) operations. Imaginary <f> 
leads to SU(2) matrices, which correspond to energy- 
preserving (passive) manipulations of the two modes or 
two states. Parametric amplifiers, four-wave mixers and 
phase-conjugating mirrors are SU(1,1) devices which can 
operate also in the quantum regime 0: they could be 
used to build a network of active (quantum) optical com- 
ponents. The corresponding passive networks could be 
realized, for example, with beam splitters or fibre cou- 
plers. Also a network of intersecting energy levels can be 
described by a network of the type in Fig.l: the avoided 
crossings between the levels are identified with the op- 
erations A and B. Corresponding physical systems are 
for instance Rydberg atoms pi and longitudinal electro- 
magnetic modes in a cavity |1. One can also consider 
the matrix P as a set of operations done by a quantum 
computer For simplicity in the following we call 

the SU(1,1) components amplifiers and the SU(2) com- 
ponents beam splitters, but actually mean any of the 
possible realizations. 

Note that according to (|l]) A — > 7 2 when the angle 
cj> — > 00, and B — > I2 when cf> — > 0. That is, in both 
of these limits the network decomposes into sets of non- 
interacting modes. Thus P describes a quantum network 
where only nearest neighbours interact, and where a sin- 
gle parameter <j> determines the relative importance of the 
interactions, i.e. the network character of the system. 

The transitions in the network are determined by the 
eigenvalues of P. The only problem in diagonalizing 



P is the relative shift between the sets of A and B. 
This can be solved by the discrete Fourier transform 
(F N ) H = exp(i2irkl/N)/VN = lu m /VN, because the 
Fourier transform of the shift matrix (Sn)m = <Jfe+i,z + 



~ 1 Ski- The whole network matrix P thus decomposes 
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and C = cosh and S = i sinh. Most textbooks present 
the solution of the Ising model in a slightly different form, 
but we have formulated the problem as in in order to 
make a connection to interferometers. The usual Mach- 
Zehnder interferometer affects the input states by a uni- 
tary transformation Um-z which can be formally written 
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where K n = exp(i<fin) is determined by a chosen phase 
<f>. Thus the Ising-type network we consider acts like 
an A^-dimensional interferometer where, instead of one- 
mode phase shifts, two-mode rotations are performed in 



between the iV-dimensional mixers F^ and Ft, |12 



From (||) one obtains the eigenvalues e ±7 ™, where j n 
are determined via 
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An explicit expression for 7„ is given via the integral 
representation H 
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As can be seen from above, the zeroth eigenvalue (n = 0) 
is not a smooth function of at 70 = 0, i.e. when 
cosh</> = y/2, its derivative ^ has a discontinuity. In 
the Ising model this gives the transition temperature 
kT c = 2.269e. We are now at the point to interpret 
what this mathematical behaviour means physically in 
the case of quantum networks. We will first consider ac- 
tive SU(1,1) networks, then the passive SU(2) ones. 
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For the active SU(1,1) networks the amplification G of 
the single components A and B is cosh 2 6 and cosh 2 <f>, 
respectively. The critical amplification G c = cosh 2 8 C = 
cosh 2 <p c = 2 has an interesting physical interpretation. It 
has been shown jl3| that for parametric amplifiers G = 2 
sets a borderline between quantum and classical perfor- 
mance of the device. For G > 2 an initially squeezed 
input loses the squeezing, i.e. its non-classical proper- 
ties in the process of amplification. This is sometimes 
called the "magic cloning limit"; if it did not exist, one 
could reproduce quantum states, which would simply vi- 
olate quantum mechanics. By considering what happens 
at the critical amplification point, as well as below and 
above it, we can show that G c coinciding with the "magic 
cloning limit" is not a mere coincidence. 

In the Ising model the transition point divides regimes 
of order and disorder. To see whether G c imposes any 
such boundary we consider again the matrices K n , now 
written in the form 
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In the limit 4> — > the matrices K n are clearly func- 
tions of u> n , but when (j) — > oo, they become increas- 
ingly independent of n; actually G < 2 corresponds 
to the condition \u) n \ > \smh(p\. Considering Eq.(j^) 
one sees that when the input is any vector (specified 
by n) of the Fourier transform matrix, i.e. of the type 
l/VN[l,l,u n ,u n ,u 2n > (jJ 2n ...], it will be affected by the 
corresponding K n matrix. For G >> 2, however, the K n 
tend to be independent of n; this means that the net- 
work gives the same response independent of the relative 
phases w" of the input modes. The opposite is true for 
G < 2. Thus the network behaves like a phase sensi- 
tive, i.e. a quantum device below G c — 2, and classically 
above it: it is logical that G c = 2 coincides with the 
quantum-classical border of the individual components. 

Let us now consider the passive SU(2) networks 
The analogy to the Ising model is then not one-to-one. 
For example ([l]) is true only for one trivial choice for the, 
now imaginary, angle <f>. We can, however, define a net- 
work which has the same basic properties as the active 
ones: nodes connected by nearest neighbour interactions, 
with one parameter <f> quantifying the importance of these 
interactions. Let us, for example, fix 6 — in/ 4 and de- 
note (j)' = —i(j>. In the limits </>'—>() and <fr' — > tt/2 
the network decomposes into sets of independent modes, 
while intermediate values of <f>' describe a network of in- 
teracting modes. The zeroth eigenvalue of the matrix P 
is defined by the equation cos 70 = ^=(cos</>' + sin0'). It 
reaches the value zero for <j)' c = 7r/4, and its derivative 
with respect to <p' has a singularity at this point. The 
transition thus takes place at the point when all the beam 
splitters are half-transmitting. Here it is, however, not a 
transition between quantum and classical regimes like in 



the case of active networks. Names for the two regimes 
separated by 4>' c can be given by considering a network of 
avoided crossings. Imagine that in Fig.l the individual 
elements are actually avoided crossings between inter- 
secting energy levels. At each crossing the system can 
either follow the energy level adiabatically, or make a 
Landau-Zener transition to the neighbouring level, that 
is, to show diabatic behaviour. Indicating which of these 
processes is more likely, we call the regime <fi' < it /4 dia- 
batic and <p' > 7r/4 adiabatic. Thus we have shown that 
the system does not evolve smoothly from the adiabatic 
to the diabatic regime and vice versa, but exhibits at the 
critical transmittance t c = \ a transition which is asso- 
ciated with singularities in the global parameters of the 
network. 



We have now identified the regimes of behaviour sepa- 
rated by the transition, both for the SU(1,1) and SU(2) 
networks. The essential feature is the appearance of sin- 
gularities at the transition point. In a normal situation 
the tuning of local parameters, that is the parameters 
of the network nodes such as beam splitters or ampli- 
fiers, leads to a smooth change in the global properties 
of the network. At the transition point this is not true. 
To demonstrate the observability of this phenomenon we 
now consider the response of the network in the case of 
two generic types of input states: the equal superposition 
of all modes, and the eigenstate of one mode. 

The output of the network is given by the transfor- 
mation P M . The Mth powers of K n in Eq.(||) can be 
obtained using the diagonalized form of (||): (K^f)u = 
coshM7„ + i sin (2irn/N) sinh M^ n / sinh 7„ = (2^)22 
and {K*f) n = (K n ) 21 sinh M 7n / sinh 7 „ = {K*f)\ 2 . An 
input in an equal superposition state with a phase peri- 
odicity determined by n, that is, [1, 0, uj n , 0, uj 2n , ...]/y/N, 
will be transformed into the output state 
\{K™) u , {Kt?) u , u" (K™ ) n , , ...]/yN. By 

choosing n one can thus control constructive and de- 
structive interference in the network (7 n +i > In)- The 
choice n = is a special one: at G c 70 = and the 
whole network becomes transparent. The transparency 
remains true independent of M; one may consider this 
to be analogous to the appearance of long range correla- 
tions at T c in the Ising model. Furthermore, for n = 
the transition is clearly manifested in the measurable 
output intensities. The global amplification coefficient 
is then simply M70, and its rate of change with respect 
to the local amplification coefficient <j>, M^r, has a 
discontinuity at G c . For an input in the eigenstate of 
the 2n + 1th mode, i.e. [0, 0, 1, 0], the 2n + 1th 
output amplitude has the form i X)m=o coshM7 m + 
i sin(27rm/W) sinh M~f m / sinh 7 m . Since this sum con- 
tains all 7 m , the singularity is smoothed. By choosing 
the type of the input one can thus modify the manifes- 
tations of the transition in the output. Similar consider- 
ations can be carried out to find the fingerprints of the 



3 



transition in the case of SU(2) networks. 

Note that due to the non-smoothness of 70 operations 
defined by G c or t c may be unstable ones. This has inter- 
esting consequences: for instance, some operations to be 
performed by the proposed quantum computer may hit 
an unstable transition point of the whole system. Also 
quantum correlations and entanglement are expected to 
show special behaviour near the transition. Furthermore, 
it could be illuminating to consider higher dimensional 
quantum interference experiments in connection with the 
Ising model. The theoretical description of networks of 
avoided crossings 15 |^] resembles the one presented 
here, with the choice A = B and with the addition of 
extra phase factors due to free evolution of the system 
between the avoided crossings. However, when the free 
evolution phase factors are multiples of n/N, the system 
reduces to the one described here, up to trivial differ- 
ences. It is interesting to note, that experimentally ob- 
served recurrence phenomena appear in these networks 
exactly for such values of the free evolution phase, that 
is, when the system reduces to an Ising-type network. 
Moreover, approximate analytical results predicting the 
recurrences are possible to derive in the adiabatic and 
diabatic limits 0], but not for the case of 50:50 crossings 
which according to our results is the transition point be- 
tween these regimes. 

In summary, we have taken a novel point of view to the 
networks used in quantum computation, communication 
and interference experiments: we have shown that they 
can exhibit transition phenomena. Although the indi- 
vidual quantum components at the network nodes are 
smoothly behaving, in certain network configurations — 
like the 2-D Ising model analogue considered here — they 
show non-smooth global behaviour. For networks based 
on active SU(1,1) components the transition is particu- 
larly interesting since the critical amplification G c = 2 
coincides with the 'magic cloning limit' above which an 
initially squeezed input loses its non-classical properties 
in the process of amplification. We have shown that the 
whole network too has regimes of quantum and classical 
behaviour: below the transition point G c = 2 the net- 
work is sensitive to the phase information in the input, 
for higher values of G the output depends mainly on the 
intensities of the inputs. In the case of passive SU(2) net- 
works the critical transmittance t c = i separates regions 
of adiabatic and diabatic behaviour. We have indicated 
how the transition is reflected in observable quantities. 
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